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We study junctions of single-channel spinless Luttinger liquids using bosonisation. We generalize 
earlier studies by allowing the junction to be superconducting and find new charge non-conserving 
low energy fixed points. We establish the existence oi g ^^ 1/g duality (where g is the Luttinger 
Liquid parameter) between the charge conserving (normal) junction and the charge non-conserving 
(superconducting) junction by evaluating and comparing the scaling dimensions of various operators 
around the fixed points in both the normal and superconducting sectors of the theory. For the 
most general two-wire junction, we show that there are two conformally invariant one-parameter 
families of fixed points which are also connected by a duality transformation. We also show that 
the stable fixed point for the two-wire superconducting junction corresponds to the situation where 
the crossed Andreev reflection (an incoming electron is transmitted as an outgoing hole) is perfect 
between the wires. For the three-wire junction, we study, in particular, the superconducting analogs 
of the chiral. Dp and the disconnected fixed points obtained earlier in the literature in the context 
of charge conserving three-wire junctions. We show that these fixed points can be stabilized for 
g < 1 (repulsive electrons) within the superconducting sector of the theory which makes them 
experimentally relevant. 

PACS numbers: 71.10.Pm,73.21.Hb,74.45.+c 
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I. INTRODUCTION 

Recently, Y-junctions of several quasi one-dimensional 
(1-D) quantum wires (QW) have been realized experi- 
mentally in single- walled carbon nanotubes [l|, I4I . Junc- 
tions of this kind are of importance for potential appli- 
cation in the fabrication of quantum circuitry. Theoret- 
ically, junctions of QW have been studied from several 
points of view 

El, 113, M M 

isation, weak interaction renormalisation group (WIRG) 
methods, conformal field theory and functional renormal- 
isation group methods. The junction has also been var- 
iously taken to be enclosing a flux, having a resonant 
level, having a Kondo spin and having a superconductor 
using one or the other techniques mentioned above. 

A comprehensive study of the junctions of three 
QW enclosing magnetic flux was carried out by Cha- 
mon et. al [y, |7|, where the wires were modeled as sin- 
gle channel spinless Luttinger liquids (LL) and confor- 
mally invariant charge conserving boundary conditions 
were identified in terms of boundary bosonic fields which 
had correspondence with a host of fixed points in the 
theory. However, superconducting junctions of multiple 
1-D QW have not been studied in the past for the case 
of arbitrarily strong electron-electron interactions. 

In this article, we study transport across multiple wires 
connected to a superconductor as depicted in Fig. [1] 
In the sub-gap region, normal refiection and transmis- 
sion of the electrons cannot occur, since charges can en- 
ter and exit the superconductor only as a Cooper pair. 
But due to the proximity effect, two new processes can 



occur. One is the phenomenon of Andreev reflection 
(AR) in which an electron like quasi-particle incident on 
normal— superconductor (NS) junction is reflected back 
as a hole along with the transfer of two electrons into the 
superconductor as a Cooper pair. The second even more 
interesting process is 'crossed Andreev reflection (CAR)' 
271 l28l |29| . |30| , whereby an electron from one wire pairs 
with an electron from another wire to form a Cooper pair 
and jumps into the superconductor, emitting a hole in the 
second wire (note that for a singlet superconductor, the 
two electrons have opposite spins). This can take place 
provided that the distance between the two wires L is 
less than or equal to the phase coherence length of the 
superconductor. Thus, for an incident electron, holes are 
either reflected or transmitted across the junction, and 
total current conservation is taken care of by the Cooper 
pairs jumping into the superconductor. However, as far 
as the multiple wire system is concerned, current is not 
conserved. The system is modeled as several 1-D LL con- 
nected to a superconducting junction. We assume that 
the width of the superconductor between any two wires 
L > a, where a is the phase coherence length of the super- 
conductor. For simplicity, we assume that the supercon- 
ductor is a singlet. Thus spin is conserved in transport 
across the superconductor and we can confine our study 
to spinless LL. For this system, we see that the super- 
conductor can be modeled simply as a (charge- violating) 
boundary condition on the bosons in the wire. We also 
find a rich fixed point structure that generalizes the ear- 
lier structure of fixed points found when multiple wires 
are connected to a normal junction. 

The superconductor explicitly violates charge conser- 
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FIG. 1: Multiple wires connected to a superconducting junction within the phase coherence length a {L < a) of the supercon- 
ductor. The processes corresponding to an incident electron on one wire undergoing CAR (amplitude represented as t^) and 
direct transmission (amplitude represented as t) to a difTerent wire across the superconductor are depicted in the figure. 



vation at the boundary, thereby it allows for a gener- 
alization of the study of Chamon et. al. to the charge 
non-conserving sector. We find that there exists a "nor- 
mal junction— superconducting junction (N— S)" duality 
given by g ^^ \/g {g is the LL parameter) between 
the charge conserving (normal) and the charge non- 
conserving (superconducting) sectors of the theory for 
junctions of any number of QW. As a consequence of this 
duality, many of the fixed points that were unstable for 
the normal junction for 5 < 1, turn out to have stable su- 
perconducting analogs. The stability of the fixed points 
mentioned here are calculated with respect to perturba- 
tions which are within the normal sector if the fixed point 
is in the normal sector and within the superconducting 
sector for the fixed point in the superconducting sector. 
The main results obtained in this article in the context 
of two-wire and three-wire junctions are : 

(a) For the most general two-wire junction, we show 
that there are two conformally invariant one- 
parameter families of fixed points which are con- 
nected to one another via a duality transformation. 
We also show that the normal sector and the super- 
conducting sector of the theory correspond to two 
distinct points on each of two one-parameter fami- 
lies of fixed points. Hence other than these special 
points on the two one-parameter families, loosely 
speaking, the rest of the fixed points represent semi- 
normal (semi-superconducting) junction. We find 
that the stable fixed point within the superconduct- 
ing sector of the theory corresponds to a situation 
where an incoming electron is completely transmit- 
ted as an outgoing hole, as shown in Fig.[2I^b). This 
is the crossed Andreev reflection (CAR). This fixed 
point is shown to be dual to the unstable connected 
(perfectly transmitting) fixed point of a two-wire 



normal junction due to the N— S duality. 
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FIG. 2: Stable fixed points of (a) normal junction (electron is 
completely reflected) and (b) superconducting junction (elec- 
tron is perfectly transmitted as a hole). 

(b) For the three-wire junction, we restrict our study 
to the special cases of normal and superconducting 
sectors. Within each sector, the theory of the three- 
wire junction effectively reduces to the most general 
theory of the two-wire junction as in both cases 
it is a theory of two independent bosonic fields. 
For the three-wire junction, out of the three in- 
dependent bosonic fields, one is pinned either by 
the charge conserving (normal) boundary condi- 
tion or by the charge non-conserving (supercon- 
ducting) boundary condition leaving behind only 
two independent fields. Hence for the three- wire 
superconducting junction also, one gets two con- 
formally invariant one-parameter families of fixed 
points, which are connected to one another via a 
duality transformation. Of all these fixed points for 
the system of a superconducting three-wire junc- 
tion, we shall mainly focus on two which are of 
interest both theoretically and experimentally : 

(i) SDp fixed point : This fixed point represents a 
junction with Z3 symmetry between the three- 
wires having maximal CAR between any of the 
two- wires. In other words, this is a fixed point 



where an incoming electron has non-zero com- 
ponents on all three wires as outgoing states. 
—2/3 of the charge is transmitted on the two 
other wires (hole transmission) and 1/3 of the 
charge is back-scattered (electron reflection). 
Note that the net change in charge at the 
boundary is e - (-2/3 - 2/3)e + (l/3)e = 
2e. This can be identified as the charge non- 
conserving analog of the Dp fixed point found 
in Ref. lZI. The SDp fixed point is shown to be 
stable for g < 1/3 within the superconducting 
sector and is identified as dual of the charge 
conserving Dp fixed point via the N— S duality. 
These fixed points are shown in Fig. [3l 





FIG. 3: (a) Dp fixed point for the normal junction. 2/3 charge 
is transmitted on each of the other wires and —1/3 charge is 
reflected, and (b) SDp fixed point for the superconducting 
junction. —2/3 charge is transmitted on each of the other 
wires and 1/3 charge is reflected. Note that we have consid- 
ered incoming electrons only along one of the wires. 



(ii) Sx± fixed point : These two fixed points, 
Sx+ and Sx-, represent a superconducting 
three-wire junction with maximally asymmet- 
ric inter-wire CAR with broken time reversal 
symmetry. An incoming electron along wire 1 
is transmitted as a hole in wire 2, and so on, 
cyclically, as shown in Fig. [3)Jb) (or the other 
way around). They are the superconducting 
analogs of the chiral fixed points (Fig. [3)^ a)) 
found earlier [a, [3, |9(] . Unlike their charge con- 
serving analogs, these fixed points are stable 
for 1/3 < g < 1. As can be seen from the sta- 
bility window of Sx± , these are the most rele- 
vant fixed points from the experimental point 
of view as they can be stabilized even for a 
very weakly interacting (g ^ 1) electron gas 
provided the charge conserving perturbations 
are weak enough. 

An extensive study of the renormalisation group evo- 
lution of several wires connected to a superconductor 
was carried out very recently in Refs. [lOl and [llj by the 
present authors and A. Saha, where conductances were 
studied in the Landauer-Buttiker language of transmis- 
sion and reflection of electrons. Interactions were taken 
into account perturbatively using the weak interaction 




FIG. 4: (a) x± and (b) Sx± flxed points. 



renormalisation group WIRG method. But for arbitrar- 
ily strong inter-electron interactions, one needs to use 
bosonization. Also, since the WIRG procedure is essen- 
tially a one-particle approach, it could only access those 
fixed points that could be expressed linearly in terms of 
fermions. To access other fixed points, one needs to use 
the technique of bosonisation. Note that some of the 
fixed points obtained from the fermionic WIRG method 
can be identified with some of the fixed points obtained 
using bosonization by taking the g close to unity limit 
but in general this is not true. In this paper, our aim 
is a comprehensive study of the system of multiple wires 
connected to a superconductor, and to identify the vari- 
ous fixed points of the system, only some of which were 
obtained in the earlier approach. 

In what follows, we first describe bosonization of su- 
perconducting junction of N number of LL wires in Sec- 
tion II. In Subsections II A, II B and II C, we apply this 
method to single-wire, two-wire and three-wire junctions 
and calculate scaling dimensions of various operators in 
the theory. In Subsection II D, we give an expression 
for conductance and calculate it for various fixed points 
obtained in previous subsections. Finally, we conclude 
with a discussion on general issues related to the physics 
of LL junctions in Section III. 



II. BOSONIZATION OF THE 
SUPERCONDUCTING JUNCTION OF LL QW 

The (spinless) electron field can be written in terms of 
bosonic fields as. 



tp(x) — Fne^'^"'" e^^^^^^^"^^^^^ + Fre"^'^'''" e^^'^^^'' 



-e{x)] 



where Fq and Fj are Klein factors for the outgoing and 
incoming fields respectively, and (/"(a;) and 9{x) are the 
dual bosonic fields and kp is the Fermi momentum. The 
wires are modeled as spinless LL on a half-line {x > 0) 
i.e., here we use a folded basis for describing the junction 
such that all the wires lie between x = and a; = oo and 
the junction is positioned at x = 0. Hence the action is 
given by 
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where prime (dot) stands for spatial (time) derivative 
and <pi{x,t) = {(pio + 4'ii)/'^, di{x,t) = {(j)io - (t'ii)/'^; 
4>i = (^/s) ^'i s-nd 9i — (vg) (f>[. 4>io and (pu are the chiral 
outgoing and incoming bosonic fields. 

The action can also be written in terms of purely the 
(fii fields or the 9i fields and as is well-known, the two 
actions are identical with the replacement of g *-> l/.g. 
The total densities and the currents in each wire can also 
be written in terms of the incoming and outgoing fields : 
the density p — po + Pi with po/i = ±(l/27r) ipQij and 
the current j = jo-ji with jo/i = ±u_F(l/27r) (f)'^^^. To 
complete the theory, the action needs to be augmented 
by a boundary condition at the origin which represents 
the physics at the junction. 

Now following the method we used in Ref. M, it is possi- 
ble to represent the junction in terms of a splitting matrix 
i.e., we connect the incoming and the outgoing current 
fields as jio\x=o = 'Cik jki\x=o |3J through a current 
splitting matrix, C^fe. For charge conserving fixed points, 
the net current flowing into the junction must be zero. 
Hence all charge conserving fixed points must satisfy the 
constraint that J2iJi — 0- ^^ terms of the bosonic fields. 

Now we can write 



this implies that ^^ 



E. 



''iO- 



a field splitting boundary condition as 

0jO = Cij(j)jl 



(2) 



which is consistent with jio — Cik jki and ^^ (p'^j — 
— J2i 4''iO ■ While writing the field splitting relations in 
the above form, we have ignored possible integration con- 
stants, inclusion of which makes no difference to the eval- 
uation of scaling dimensions of operators around various 
fixed points. Hence, the current splitting and the field 
splitting matrices are taken to be the same. 

For the boundary condition or equivalently for the ma- 
trix Cy to represent a fixed point, it should not flow un- 
der RG. This means that it has to be scale invariant or 
equivalently in 1 -I- 1 dimensions, conformally invariant. 
Here, this condition simply means that the trace of the 
energy-momentum tensor of the bosonic fields has to van- 
ish at the boundary (x = 0). This yields the following 
condition Isil |32| 



^^i f[\x=0 = 



(3) 



where, the (pi are mutually independent fields such that 
for a junction of A^-wire system with N' constraints. 



771 — N — N' . Hence, for two or more fields coupled to 
the junction, one can have mixed boundary conditions, 
besides the Dirichlet (</>, = 0) and Neumann {(p^ = 0) 
boundary conditions. For m = 2 (which is the case for 
the most general two- wire junction or the three- wire case 
in either the purely charge conserving or the purely su- 
perconducting limit), there are two independent families 
of solutions possible to the above equation given by 



(4) 
(5) 

where a is a real constant, independent of x and t. In 
terms of the 2x2 field splitting matrices, it is easy to 
check that this is equivalent to taking 
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-aV2, 


and 


<P2 = O-'Pl 
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-0.1^2, 
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ip2 = aip[ 
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- Sl Cl 



and 
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C2 S2 
S2 -C2 



(6) 



where c; = cos 6i and Si = sin 9i are real parameters and 
Cgj and Cgj are the field splitting matrices at the junc- 
tion for the if fields. The two family of solutions are con- 
nected via duality transformation, which we call Dg^ej. 
DsiSj duality can be accomplished by either ipi ^^ di 
or ip2 <-^ '(?2- For the two- wire system, tpi and 'di can 
be identified with the (pi and 9i respectively, i being 
the wire index and C^^ and <Ce2 can be identified with 
the current splitting matrices for the two-wire system. 
For the three-wire case, (pi and di have to be taken 
to be linear combinations of (jj^s and 9iS after impos- 
ing normal and superconducting boundary conditions. 
For these one-parameter families it turns out that the 
incoming and outgoing (bosonic) boundary fields sat- 
isfy the bosonic commutation relations of the bulk given 
by [0o//(2;)>o//(a;')] = ±iTTSgn{x - x'), so imposing 
bosonic commutation relations gives no new constraints. 
The boundary conditions may also be written in terms 
of the Boguliobov transformed free bosonic fields, (f>io/i 
which are defined as [9| 
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{g + I) ^^o - [g - I) ^^I 
-{g- l)cf>io + {9 + l)0i7 



(7) 
(8) 



For the tilde fields, the boundary condition cpio = Cij <pji 
for the A^-wire junction becomes 



^ij (Pjl 



with 



[(g + l)C+(g-l)I] 
[(5-l)C + (.g + l)I] 



(9) 



Thus M is the matrix that connects 'free' incoming and 
outgoing bosonic fields whose dimensions we know how 



to compute. Now notice that when C^ = I, the above 
equation simpUfies to (pio = Cij 4>ji, but not otherwise. 
This impHes that, for the case of C^ =1 both the inter- 
acting fields {(pio/i) and the free fields {4>io/i) satisfy the 
same boundary condition. Also note that current conser- 
vation implies that the elements of the splitting matrix 
C are real and satisfy the constraint, 



E< 



(10) 



Furthermore the constraint that both the incoming and 
outgoing fields satisfy bosonic commutation relations in- 
dependently implies [9| 



E< 



= 1 



and 2_, ' 



■'ij^i+lj 



= 



(11) 



which is essentially the same constraint that is obtained 
from imposing the constraint of scale invariance, or re- 
quiring Cij to be a fixed point. 

For the three-wire system, most of the fixed points 
studied in Ref. [vl can be obtained as C matrices satisfy- 
ing the above constraints. For instance, the disconnected 
normal (DN3) fixed point where each of the wires inde- 
pendently has a Neumann boundary condition on the (fi 
field at origin corresponds to C = I and the Dp fixed 
point has C matrix of the form 



- 1/3 2/3 2/3 
2/3 -1/3 2/3 
2/3 2/3 -1/3 



(12) 



It turns out that several other C-matrices obeying 
the constraints mentioned above fall into the two one- 
parameter families given in Eq. [S] and hence can be iden- 
tified as conformally invariant fixed points. Also note 
that both the disconnected DN3 and the above Dp fixed 
points belong to the special class of Cij matrix for which 

C2 =1. 

Physically, the disconnected DN3 fixed point (called N 
in Ref. 7) corresponds to a situation where the conduc- 
tance between any two-wires is zero, whereas the Dp fixed 
point corresponds to a situation where there is a perfect 
symmetry among the three-wires and the conductance 
between any two- wires has the maximal value allowed by 
Z3 symmetry. Note that this maximum is larger than 
the maximal inter-wire conductance that would be al- 
lowed within a free-electron picture for the maximally 
conducting Z3 symmetric case [7[ and this is related to 
the fact that for the bosonic Z3 symmetric fixed point, 
multi-particle scattering leads to an enhancement of con- 
ductance as was discussed in Ref. l7|. In Subsection II D, 
by calculating the conductance, we will show that for the 
analogous situation in the superconducting sector, this is 



no longer true, and in fact, there is a reduction in the 
conductance as compared to the free electron case. The 
difference in the processes participating in the two sectors 
can also be seen in Fig. [3l 

The charge conserving constraint at the junction im- 
plies that the boundary condition on the CM field de- 
fined as (j)cM = '^i4>i always has to be Neumann i.e., 
S 4>ii = S 4'iO + c, where c is a constant. However, 
in the presence of a superconducting junction strongly 
coupled to the wires, there will only be charge non- 
conserving processes at the boundary (i.e., it can either 
absorb or emit a Cooper-pair), and charge conserving 
processes will be suppressed (at energies below the super- 
conducting gap). Now if we impose Dirichlet boundary 
condition on the CM field, it turns out that it gives the 
correct boundary condition at the junction that converts 
an electron to a hole and vice- versa, and mimics the exis- 
tence of a superconductor at the junction. This leads to 
new fixed points, which have not been explored in Ref. |7|. 
This is one of the main points of our article. 

In order to establish the duality between the normal 
and superconducting junctions, let us consider the case 
of an NS junction where a single QW is connected to 
a superconductor (This case was considered briefly in 
the appendix of Ref. 7) in the sub-gap regime. In the 
limit, when the coupling between the wire and the su- 
perconductor is strong (i.e., there is no back-scattering 
of electrons) , the system is in the perfect Andreev limit 
and hence an incoming electron current is completely re- 
flected as an outgoing hole current i.e., ji ~ —jo- We caU 
this as the Andreev (Ai) fixed point. This implies that 
the boundary condition on the (/)(x = 0, t) field is Dirich- 
let (or equivalently Neumann on the dual 9[x — 0,t) 
field) and the total current at the junction is given by 
j = ji — jo = 2jj . This can be easily generalized to a 
system of superconducting junction of A^-wires. 

For the A^-wire system, we must have the sum of the 
incoming electron current equal to the sum of outgoing 
hole current, which means that "^^jn + ^ijio = at 
the junction. In turn, this implies that ^^ Pi — 0, i.e., 
the total electron density is zero at the junction. This 
is of course the correct boundary condition as the elec- 
tron density is expected to vanish at the junction due 
to the finite superconducting gap. In terms of the split- 
ting matrix C, the above constraint translates into the 
condition, 



E' 



-1 



(13) 



in contrast to the current conserving constraint 
(Eq. [T0|) [35|. The other constraints coming from the 
bosonic commutation relations that (pio/i have to sat- 
isfy, given by Eq. 1111 still remain valid. As mentioned 
earlier, these matrices fall into the two one parameter 
families given in Eq. [6l thus enabling us to identify them 
as fixed points. In fact, given a Cij matrix representing 



a fixed point in the normal sector, its dual fixed point 
in the superconducting sector can be obtained by trans- 



forming Ci 



~Cij. It can be easily checked that this 



prescription of finding the dual fixed points is consistent 
with the constraints given by Eas.[TUland[T51 

The duality between the charge conserving and the su- 
perconducting boundary conditions is now obvious and 
can be understood physically as follows. Current conser- 
vation implies that the net current should be zero at the 
junction while, superconductivity implies that the net 
electron density at the junction has to be zero, due to 
the existence of the gap for single electron excitations 
in the superconductor. So, in the current conserving 
case, the boundary condition on the (j)cM held is Neu- 
mann (or Dirichlet on 9cm field, i.e., J2iJii^^^) ~ 0) 
while for the superconducting case, the boundary condi- 
tion is Dirichlet on 0cm field (or Neumann on ^cm field, 
i.e., J^iPii^T^) — 0)- As the 9 and the (f> fields have 
5 <-> l/g duality among themselves, it automatically ex- 
tends to the various fixed points in one sector and their 
analogs in the other sector, which are obtained by impos- 
ing further boundary conditions on the fields other than 
the CM field. We confirm this by explicitly calculating 
the scaling dimension of operators corresponding to all 
possible perturbations around these various fixed points. 

Note however that the N— S duality exists over and 
above the dualities that exist within each sector. For in- 
stance, within the charge conserving sector for the two- 
wire system, there exists a duality between weak back- 
scattering (strong tunneling) and strong back-scattering 
(weak tunneling) limits with g <-^- 1/g interchange. Sim- 
ilarly in the superconducting sector also, there exists a 
duality between weak back-scattering of holes or weak 
Andreev refiection (strong transmission of holes or strong 
CAR) and strong back-scattering of holes or strong An- 
dreev reflection (weak transmission of holes or weak CAR) 
with g <-> 1/g. This essentially follows from the Dg^^g^ du- 
ality. 

We will now explicitly consider the cases where there 
are N = 1,2 and 3 wires coupled to the superconductor. 



A. Single-'wire junction 



■^- 



■^- 



■^- 



(a) 



(b) 



FIG. 5: (a) Andreev reflection and (b) normal reflection from 
a superconductor. 



We start with the simplest case of the NS junction 
where the number of wires, A^ = 1. In this case, there 
are two single element splitting matrices that satisfy the 
constraints of Eq. [Til and only one of them satisfies the 
superconducting constraint of Cy = — 1. In that case the 
wire is perfectly connected to the superconductor and an 
incoming electron is scattered back perfectly into a hole 
(see Fig. [5]) . This is the perfect Andreev limit described 
before where jn — —jto- The scaling dimension of the 
electron back-scattering operator, 'ipj4'0 (the subscripts 
I/O on the electron fields refer to incoming and outgoing 
branches) around this fixed point can be easily found 
by bosonizing it as i/'/'0o ^ e^^'^^e*'^'^. Upon writing 
it in terms of the Boguliobov transformed fields, we can 
compute the scaling dimension of this operator to be 2g. 
Note that the back-scattering operator we have turned on 
around the charge non-conserving fixed point is charge 
conserving. 

The other fixed point corresponds to the charge con- 
serving case, where the splitting matrix is C^j = -1-1. 
Here the incoming current is perfectly (normal) refiected 
{jii — jio) and the wire is completely disconnected from 
the superconductor. We can now turn on a charge vio- 
lating perturbation, such as the Andreev refiection (AR) 
operator, V'/V'o ~ e^'^^e"^". The scaling dimension of 
this operator turns out to be 2/g. This establishes the 
9 ^^ ^l9 N— S duality between these two cases. 



B. Two-wire junction 

Let us now go on to case of the NSN junction, where 
the number of wires is A = 2. In this case, the cur- 
rent splitting matrix is 2 x 2. Unlike the previous case 
(NS junction), here we find that there are two fixed points 
in the superconducting sector and they are represented 
by the following two matrices 



'1 — 



-1 
-1 



and 



-2 — 



-1 
-1 



(14) 



The matrix Ci corresponds to a situation where the two- 
wires are individually tuned to the disconnected Andreev 
(A2) fixed point (electrons are reflected back as holes) 
whereas the matrix C2 implies perfect CAR between the 
wires and is called the crossed Andreev (CA2) fixed point 
(electrons perfectly transmitted as holes). As can be eas- 
ily checked, Ci is a particular case of Ce^ (6*1 = tt) and 
C2 is a particular case of C-g^ (^2 = — ""/S) given by Eq.[6l 
It is easy to see that these two cases are analogous to 
the completely reflecting (disconnected) and completely 
transmitting (fully connected) cases for the normal two- 
wire junction. 

Let us now turn on tunneling or back-scattering opera- 
tors as perturbations around these flxed points. Around 
Ci, which is fully disconnected, we switch on a CAR op- 
erator which will convert an incoming electron in one 



wire to an outgoing hole in another, given by 4'ii'^20 ^ 
The dimension of this operator can be com- 



ol-VlI p1-V20 



puted by re-expressing the operator in terms of the 
Boguhobov transformed fields. 

Since the matrix Ci is just the negative of the iden- 
tity matrix, it is trivial to see that the Boguliobov trans- 
formed fields also satisfy the same boundary conditions 
as the original fields. The scaling dimension can easily 
be computed and it turns out to be equal to g. Anal- 
ogously, around the C2 fixed point, where an electron 
injected in the first wire gets perfectly transmitted as a 
hole in the second wire, we can switch on the AR oper- 
ator, V'i/V'iO ~ g^<t>ii gi-tkio ^ Again the bosonic fields can 
be re-expressed in terms of the Boguliobov transformed 
(j) fields and since C2 = I, the tilde fields also satisfy the 
same boundary conditions. Here, we find that the An- 
dreev back-scattering operator has the dimension l/g. 
Hence, within the superconducting sector, for repulsive 
inter-electron interactions, g < 1, the CA2 fixed point is a 
stable fixed point (shown in Fig.[2l^b)) while the fully dis- 
connected A2 fixed point is unstable. This is in contrast 
to the normal charge-conserving junction of two-wires, 
where the "cut" wire (shown in Fig.[2Ua)) corresponds to 
stable fixed point for repulsive interactions. Again this 
can also be understood in terms of the g ^-s- l/g N— S du- 
ality. 

In the above analysis, we have restricted ourselves to 
either Neumann (for normal) or Dirichlct (for supercon- 
ducting) boundary conditions on the CM field and then 
analysed the system, which essentially reduces the system 
to a single boson (m — 1) problem. However, once we 
allow for arbitrary charge non-conservation, then for the 
two- wire system, both the CM field -75(</^i + (f>2) and the 
relative field -j^{4'i — 4'2) enter the picture. Hence, the 
system can no longer be reduced to a single boson theory 
as could be done when the charge conserving or the su- 
perconducting boundary condition removed the CM field 
completely from the scene. 

Hence in general for a two- wire junction, we have a gen- 
uine m — N — 2 problem. As mentioned earlier, in terms 
of these two fields, scale invariance of the boundary con- 
dition gives us two one-parameter family of fixed points 
which are consistent with bosonic commutation rules im- 
posed on the incoming and outgoing fields. The two fam- 
ilies are connected via duality transformation {(f> <-> 9) on 
either the (t>i field or, the (t>2 field, where 1 and 2 are wire 



indices. So in conclusion, the important point to note is 
that except for 6*1 = ("cut") and 6*2 = 7r/2 ("healed") 
for the normal case or, 61 = ir and 62 = —tt/2 for the 
superconducting case, these fixed points belong neither 
to the category of charge conserving fixed points nor to 
the category of superconducting fixed points. A similar 
isolated fixed point, called Andreev-Grifhths (AG) fixed 
point which allowed both superconducting and charge 
conserving transmissions and reflections was seen earlier 
in WIRG formalism in Refs.llOandllll by the authors and 
A. Saha. 



C. Three-wire junction 





FIG. 6: (a) DN3 and (b) A3 fixed points. 

Finally, let us consider the case where there is a su- 
perconducting junction of three LL wires, i.e. N = 3 and 
the current splitting matrix is 3 x 3. Here too, just as 
in the normal three- wire case 7] , we do not have a com- 
plete classification of all the fixed points of the system 
in general. However, a partial classification within the 
superconducting or the normal sector can be obtained in 
terms of the current splitting matrix which can be derived 
from the Cgj and Cg^ matrices given in Eq. [6l For the 
superconducting case, it is easy to see that we will have 
a fixed point corresponding to the situation where each 
individual wire is tuned to the Andreev fixed point (see 
Fig. [5]) at the junction. This is the disconnected Andreev 
(A3) fixed point. We will also have the SDp fixed point 
as mentioned earlier. These are the analogs of DN3 and 
Dp fixed points for the normal sector. We will focus on 
these fixed points first. The current splitting matrices 
representing the A3 and the SDp are given by 



J 



^3 

-^3 



- 1 











-1 











— ; 



and 



-3 

^SDf 



1/3 -2/3 -2/3 
-2/3 1/3 -2/3 
-2/3 -2/3 1/3 



(15) 



Let us now compute the stability of these two fixed points. Around A3, the CAR operator is given by 
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■0»/'0jO ^ e''*'^'e"^^° where j ^ i and i,j = 1,2,3. This 
is the same operator that we considered in the two-wire 
case around the A2 fixed point, and the dimension of the 
operator of course turns out to be g because essentially 
it only involves tunneling between two-wires that are dis- 
connected from each other. This is dual to the scaling 
dimension of the normal tunneling operator around the 
disconnected three-wire fixed point, which is l/g. Thus, 
it gives a simple check of the general N— S duality. 

A more non-trivial check is to consider the stability 
of the SDp fixed point. The general N— S duality implies 
that this should be dual to the usual Dp fixed point of the 
normal junction. Let us first consider the operator cor- 
responding to CAR i.e., i'j^ii'jQ ^ e"^'^e~"^^°. Since the 

matrix (C^ Y ^ Ij the Boguliobov transformed bosons 
also satisfy the same boundary condition as the original 
fields, as mentioned earlier. Hence the dimension of the 
operator can easily be computed to be l/Sg. Now con- 
sider the scaling dimension of the normal tunneling op- 
erator, ipijipjQ ~ e"^'^e~*'^-'° around the normal Dp fixed 
point. This has been computed in Ref.Qto be g/3. Thus 
the dimensions of these operators are related by g ^-» 1/g 
N— S duality. Similarly, if we consider the AR operator 
in each wire, 4'iii'io ~ e*'^'^e*'^'°, its dimension can be 
computed to be 4/3(7. This is dual to the dimension of 
the usual reflection operator, V'i/V'Jo ~ e*"^'^e~*'^'° in a 
normal junction which was earlier found to be ^g/2> [7|. 

Finally, let us consider the tunneling of the incoming 
electron in wire i to the incoming electron in wire j. Here, 
since the tunneling happens within the incoming channels 
before the electron reaches the junction, the operator is 
given by V'i/V'I/ ~ e^'^'^e"*'^'^. In other words, this is a 
charge conserving operator, unlike the two other (charge 
violating) operators for which we calculated the scaling 
dimensions. Hence the third tunneling operator that we 
consider as a perturbation around the SDp fixed point 
is the same as that for the Dp fixed point. However, 
the scaling dimension of this operator computed around 
the SDp and the Dp fixed points turn out to be differ- 
ent because the boundary condition explicitly enters the 
computation of the scaling dimensions. They turn out 
to be 1/g and g respectively, which is in accord with the 
N-S duahty. 

To sum up, we find that the scaling dimensions of 
the three classes of operators (V'i/V'jo; '^u'^io'^ '^u'Vji) 



around the SDp fixed point to be (1/3(7; 4/3(7; 1/ff)- 
These are connected by N— S duality to the three classes 

of charge conserving operators (■(/'i/V'jo; '^u'^lo'' '^u'^ji) 
around the Dp fixed point which turn out to be 

{g/3; 4^/3; g). This actually exhausts the set of all 
possible operators allowed by symmetry within the su- 
perconducting sector around the SDp fixed point. The 
important point to note is that for values of g such that 
g < 1/3, all these operators are irrelevant. So in the 
limit where the junction is tuned such that there are no 
charge conserving normal tunnelings or reflections at the 
junction, this fixed point is stable for g < 1/3. 

One can also perturbatively include the effect of charge 
conserving tunneling and reflection processes at the junc- 
tion and calculate their scaling dimensions. Such opera- 
tors are also connected by the same N— S duality trans- 
formations between the charge conserving and the su- 
perconducting sectors. The charge conserving tunneling 
operator, V'i/V'Io '^ e^'^^' e~^'^i' between two- wires across 
the superconducting junction around SDp fixed point is 
found to have a scaling dimension (2(7^ -I- 3)/3(7. This 
continues to be irrelevant for g < 1/3; hence it does not 
disturb the stability of the SDp fixed point. But the scal- 
ing dimension of the normal reflection around SDp flxed 
point turns out to be 2g/3, which is relevant for g < 1. 
Hence, we conclude that the SDp fixed point is stable 
within the superconducting sector for g < 1/3 but in- 
cluding the normal refiection operator will make it flow 
to the disconnected charge conserving fixed point. We 
can also view the SDp fixed point as a strong tunnel- 
ing limit of the CAR processes around the A3 fixed point. 
This is analogous to the strong coupling— weak coupling 
duality between the Dp and the DN3 fixed points for the 
normal three- wire system as was pointed out in Ref. |7|. If 
we now compare the scaling dimensions of the CAR oper- 
ator between the two flxed points, these are g and 1/3(7 
respectively for the A3 and the SDp flxed points. In con- 
trast, for the disconnected fixed point and Dp fixed point 
of the normal wire, they are 1/g and g/3 respectively. 
Hence the duality 1/g ^^ g/3 for the normal case goes 
over to (7 <-> l/3.g in superconducting case in agreement 
with the N— S duality. 

Next we consider the Sx± fixed point. This fixed point 
is described by the following two current splitting matri- 
ces given by 



^3 




and Cc 




(16) 



Here the subscript S stands for superconducting case and the +/— for chirality. The C^s+ fixed point cor- 



Normal junction 


Superconducting junction 


(a) N=l 


(a) N=l 




Disconnected fixed point DNi, stable for g < 2 


Andreev fixed point Ai, stable for g < 1/2 




(b) N=2 


(b) N=2 




Disconnected fixed point DN2, stable for g < 1 


Andreev fixed point A2, stable for g > 1 




Fully connected fixed point, stable for g > 1 


Crossed Andreev fixed point CA2, stable for g < 1 




(c) N=3 


(c) N=3 




Disconnected fixed point DN3, stable for g < 1 


Andreev fixed point A3, stable for g > 1 




Chiral fixed point, x±, stable for 1 < gr < 3 


Supercond. chiral fixed point, Sx±, stable for 1/3 < g < 


1 


Dp fixed point, stable for <? > 3 


SDp fixed point, stable for g < 1/3 





TABLE I: Stability of some of the fixed points for different values of the interaction parameter g for A'' = 1,2 and 3 wires 
connected to a normal junction and a superconducting junction are tabulated. 



responds to a situation where there is perfect CAR of 
electron from wire 1 ^ 2,2 ^ 3, 3 ^ 1. On the other 
hand, 'C^g_ corresponds to a situation where there is per- 
fect CAR of electron from wire 1 ^- 3,3 ^ 2, 2 ^- 1. 
The fixed point for one of these cases along with the 
analogous fixed point for the normal junction is shown 
in Fig. m Both these fixed points break time reversal 
symmetry and depend on the direction of the effective 
magnetic field through the junction. 

Next we calculate the scaling dimensions of all possi- 
ble operators around these fixed points, which are the 



following : (i) AR in each wire, 
(ii) CAR between any two wires, 






ai<t>il pi4>jO ■ 



and (iii) normal tunneling between the incoming chiral 
branches of any two wires, tpn^'li ~ e*'^"e~"^^^. Analo- 
gous to the normal chiral case, the scaling dimensions of 
all these operators turn out to the same and are given by 
45/(3g^ + l) for both C|+ and C|_. Note that all the op- 
erators listed above are marginal for g = 1 and g — 1/3. 
Now let us compare these with the scaling dimension of 
all possible operators around the chiral fixed point for 
the normal junction. The current splitting matrix for the 
normal case just requires the replacement of —1 by 1 for 
both the C|._^ and C|_ in Eq.[Tni The scaling dimensions 
of all possible operators that can be switched on around 
either of the fixed points are given by 4(7/(3 + g^). ft 
is easy to check that the scaling dimension of the opera- 
tors around the Sx± fixed point (represented by C^ , and 
C5_) is related to that of the operators around the nor- 
mal chiral fixed point \yy g ^^ 1/g. This is as expected 
from the N— S duality relation between the superconduct- 
ing and the normal sectors. However, unlike the normal 
chiral fixed point which is stable for 1 < g < 3 (attractive 
electrons, hence unphysical), the Sx± fixed point is stable 
for 1/3 < g < 1 (repulsive electrons); this fact makes this 
fixed point experimentally relevant as this fixed point can 
be stabilized even for very weakly interacting electrons. 

Now we consider the influence of charge conserv- 
ing operators corresponding to tunneling of electron 



tion between any two- wires and normal reflection of elec- 
trons (V'i/V'io ~ e«0i/e-«0io) within each wire. The 
scaling dimensions of both these operators are given by 
(25(1 -I- g^))/{l + 3(7^) and hence these operators are rel- 
evant for g < 1. This means that these fixed points 
are stable only within the superconducting sector for 
1/3 < g < 1 but not in general. Hence the Sx± fixed 
point can be relevant for experiments if one has weakly 
interacting electrons but for strongly interacting elec- 
trons (g < 1/3) the relevant fixed point would be the 
SDp fixed point as long as normal reflection at the junc- 
tion is reasonably weak. 

We summarize the results of this section in Table HI 



D. The conductance matrix and N S duality 

As shown explicitly in the appendix of Ref. 26, the 
conductance tensor relating the current and voltage in 
different wires can be computed from the Kubo formula 
and can be related to correlation functions of the incom- 
ing and outgoing currents. These correlation functions 
can be computed in terms of the R^ (5) matrix relating 
the incoming and outgoing free fields. Note from Eq. ^ 
that this matrix is in general a function of g. 

So far, we have discussed the dualities existing between 
the fixed points in the normal and superconducting sec- 
tors in terms of scaling dimensions of various perturba- 
tions around the fixed points. Now we can quantify the 
duality in terms of the fixed point conductances of the 
various fixed points in the normal sector and the corre- 
sponding dual fixed points in the superconducting sector. 

Let M^ represent the M matrix (which connects the 
'free' incoming and outgoing fields and which has been 
defined in Eq. [5]) for a fixed point in the normal sector. 
Then the R matrix corresponding to the dual fixed point 
in the superconducting sector is given by 



M^(5) 



pJV 



(1/3) 



(17) 



ii'iii'jo ~ e^'^^' e ^'t'io) across the superconducting June- This can explicitly be checked from Eq. [5]when we make 
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Cij . In terms of the 



the duahty transformation, C^ - 

the R^ and R"^ matrices, the conductance matrix is given 

bylH 



^N 



h 



:,N\ 



and 



56 



(I-R^)(18) 



where N and S stand for normal and superconducting 
cases respectively and the identity factor just represents 
the incoming current along each wire. For the two-wire 
(NSN) junction the stable fixed point within the super- 
conducting sector for g < 1 was found to be the A2 fixed 
point. The R matrix is identical to the C matrix for this 
case as C^ = I and hence 



Gl = ^ 




(19) 



This implies that if the two wires (labelled 1 and 2) are 
biased with respect to the superconductor at voltages, Vi 
and V2 then for an injected electron current, (e^/ft.) Vi the 
jmiction will suck in an electron current equal to (e^//i)V2 
from wire 2. Next let us consider the SDp fixed point for 
the three-wire case. Here also as C^ = I, the R matrix is 
identical to the C matrix and hence the conductance can 
be directly written down using Eq. [TSJ 






ge^ 
h 




(20) 



This is to be contrasted with the conductance matrix for 
the Z3 symmetric AG fixed point [10| for the three-wire 
system within the free-electron manifold obtained using 
WIRG method given by 



^free— electron 



\ (■ + ^) 

^2 / 10/9 4/9 4/9 

4/9 10/9 4/9 

^^ 4/9 4/9 10/9 



h 



(21) 



where, the elements of S are obtained from the S-matrix 
for the free electron problem by taking squares of the 
absolute values of individual elements themselves. This 
implies that § is just the current splitting matrix for the 
free electron case. 

In Ref. [tI, it was pointed out that for the normal junc- 
tion, the diagonal conductance for the Dp fixed point 
with Fermi liquid leads, (Gf, )ii = Ae^/3h was larger 
than its free electron counterpart 8e^/9h for the Grif- 
fiths fixed point [5(. In contrast, we find that this is 
no longer true for the superconducting case. The free 



electron conductance for the AG fixed point is lOe^/9/i, 
which is actually larger than the diagonal conductance 
for the SDp fixed point with Fermi liquid leads given by 
{^SDp)ii = e^//i. So in conclusion, electron-electron in- 
teractions lead to either enhancement or suppression of 
conductance with respect to its free electron counterpart, 
depending on whether the junction is normal or super- 
conducting. 

Finally, we consider the conductance matrix for the 
Sx± fixed point. For this case also, one can compute the 
conductance by evaluating the R matrices corresponding 
to the C|^ and C|_(Eq.[16l) which are. 






^C% 



-1+9' 


29(1+9) 


2(-l+9)9 


l+3g2 
2(-i+g)g 


1+39- 
-1+9' 


l+3g2 
29(1+9) 


l+3g2 
29(1+3) 


1+392 
2(-l+9)9 


l+3g2 

-1+9' 


l+3g2 


1+392 


1+392 


-1+g^ 


2(-i+g)g 


2g(i+g) 


l+3g2 

2g(i+g) 


1+3^2 
-l+g2 


1+392 
2(-i+g)g 


l+3g2 
2(-i+g)g 


l+3g2 

2g(i+g) 


1+3^2 
-i+g^ 


l+3g^ 


1+392 


l+3g2 



(22) 



(23) 



Now the conductance matrix immediately follows from 
Eq. [THl Unlike the case of SDp fixed point, the R ma- 
trix in this case depends on the LL parameter g. Given 
this, one can now explicitly check the duality as defined 
in Eq. 1171 between the chiral fixed points for the nor- 
mal junction and the Sx± fixed points in the supercon- 
ducting sector, by comparing the conductance matrix for 
Sx± fixed points with the conductance matrix obtained 
in Ref. 'T" for the normal chiral fixed points. 

Hence with all the explicit checks via the scaling di- 
mensions of operators and conductance calculations, we 
have established N— S duality rigorously for junctions 
comprising of one, two and three-wires. In general, the 
N— S duality holds for junctions of any number of wires 
including those with A^ > 3. 



III. DISCUSSION 

In the context of three- wires, both within the supercon- 
ducting sector and in the normal sector, there are only 
two independent bosonic fields connected to the junc- 
tion, since the CM field is constrained either by charge 
conservation or by the superconducting boundary condi- 
tion. As in the superconducting two-wire junction, this 
implies that scale invariance of the boundary condition 
leads to two one-parameter families of fixed points, which 
are related by a Dg^gj duality transformation on one of 
the two independent fields. But for the three wire system 
the Dgj^g^ duality transformation is not a one to one map 
between the two one parameter families of fixed points. 
For example, for the superconducting three-wire junc- 
tion, when the Dg^^g^ duality transformation is made on 
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the -75(01 — 02 ) field, the A3 fixed point (a representa- 
tive of the C^i family) goes to the SDa fixed point (a Z3 
asymmetric fixed point where two of the three wires are 
stuck at the CA fixed point and the third wire is tuned to 
the Ai fixed point with the junction) fixed point which is 
in the Cg^ family. But when the duality transformation 
is made on the -7s(0i + 02 " 2^3) field, the disconnected 
fixed point goes to the fixed point given by the matrix 
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^3 



-2/3 1/3 -2/3 
-2/3 -2/3 1/3 
1/3 -2/3 -2/3 



t 



in the Cg^ family. 

One of the most intriguing observations in the context 
of junctions of LL in general is that in the 5 — > 1 limit, 
all the fixed points in the bosonic language which have 
entries zero and unity in the current splitting matrix re- 
duce to the fermionic fixed points which have entries zero 
and unity in the current splitting matrix (obtained from 
the free electron S-matrix) which were obtained using 
WIRG formalism. The scaling dimensions of all possible 
perturbations around these fixed points also match in this 
limit. But this is not generally true when the entries are 
not and 1. For instance, the Dp fixed point in bosonic 
language is not the same as the fermionic Griffiths fixed 
point in the WIRG formalism, as neither the conductance 
at the fixed point, nor the scaling dimensions of the op- 
erators around the fixed point match in the 5 ^ 1 limit 
even though they share the same Z^ symmetry and the 
current splitting matrix for the Dp fixed point is identical 
to the §- matrix for the fermionic Griffiths fixed point. 

Also from the N— S duality, we should expect the 
analogs of the Dm and the M fixed points of Refs. [y and 
lZI to exist in the superconducting three- wire system, but 
a more detailed analysis of these fixed points is beyond 
the scope of the present work. 

Last but not the least it is very interesting to note 
that all the fixed points that we have considered in this 
article are noiseless (because there is no probabilistic par- 
titioning of the current) even though many of them do 
not correspond to the perfectly transmitting or perfectly 
reflecting situations. This is unlike the fermionic fixed 
point obtained using WIRG formalism, where the Grif- 
fiths fixed point and the Andreev-Griffiths fixed points 
had probabilistic partitioning of the current and were 
noisy. 
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